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Abstract 

The nucleus 44 Ti has low-lying levels of various kinds of mutually very different nu- 
clear structure displaying the richness of the nuclear many-body dynamics. It is shown 
| : that the deformed-basis antisymmetric molecular dynamics by the use of the Goguy 

D1S force reproduces successfully and unifiedly two types of coexistence phenomena in 
44 Ti. Namely, on one hand, the coexistence of the mean field structure and the cluster 
structure is confirmed by verifying the normally deformed structure of the K 7T =3 : [ band 
with a 1-particle-jump intrinsic configuration and the a+ 40 Ca cluster structure of the 
K^—O^ band. The mixed character of the mean-field-like structure and the a+ 40 Ca 
cluster structure of the ground band is also shown. On the other hand, the coexistence 
of the normal deformed mean field and the superdeformed mean field is confirmed 
by verifying the triaxial superdeformation of the K % =Qi^ band and the K n =2f band 
which has a 4-particle-jump intrinsic configuration. Good reproduction of the experi- 
mental data is shown for many kinds of quantities including the energy spectra, electric 
transition rates, alpha spectroscopic factors. Preliminary discussions are given on the 
existence of hyperdeformed excited states, the relation between superdeformation and 
clustering and so on. 
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§1. Introduction 



The nucleus 44 Ti is the lightest N=Z even-even nucleus in the pf shell region and its 
structure is very important for the study of the neutron-proton correlation in the pf shell 
region. Since we can now expect the increasing experimental information about the proton- 
rich nuclei along the N~Z line up to A~100, the study of the neutron-proton correlation in 
44 Ti is even more important than before. 44 Ti is analogous to 20 Ne in the sense that both 
nuclei have two neutrons and two protons outside the doubly closed-shell core nuclei 40 Ca 
and 16 0, respectively. It is well known 1 )' 2 ) that 20 Ne has non-small a-clustering character 
even in the ground band states and much prominent a+ 16 cluster structure in some of its 
excited bands. Therefore there have been many microscopic theoretical works to investigate 
possible a-cluster structure in 44 Ti. 3 )' 4 ) However, more convincing argument 5 -* than these 
microscopic studies about the a-clustering in 44 Ti was obtained from the knowledge of the 
so-called unique a- 40 Ca optical potential which was established by the study of the elastic 
a- 40 Ca scattering including ALAS (anomalous large angle scattering) and nuclear rainbow 
phenomena. 6 ^ From the fact that the band head + state of the lowest rotational band of 
this unique a- 40 Ca optical potential is located several MeV below the a- 40 Ca threshold, the 
lowest band was assigned to correspond to the ground rotational band of 44 Ti, 5 ) which means 
that the ground band states of 44 Ti should have sizable amount of component of a+ 40 Ca 
cluster structure. Furthermore, the K n —0~ band which is the parity-doublet partner of 
the lowest band ( K w = + ) was predicted to exist with its band head 1~ state a few MeV 
higher above the a+ 40 Ca threshold. A similar argument was also presented by a microscopic 
approach 7 ) which fitted a- 40 Ca scattering data in addition to the usual structure study of 
the bound states of 44 Ti. The parity-doublet partner band with K n —0~ was actually found 
later experimentally by using a-transfer reaction. 8 ) 

It is usually considered that the formation of the mean field becomes much more firm 
in the pf shell and heavier mass region compared with the smaller mass region where the 
a clustering is active especially below the mass A~30. However, recent experiments show 
that the mean field around the mass A~40 is not simply static but dynamic in its nature. It 
is because recently many low-lying rotational spectra associated with superdefomed bands 
have been observed in rather light mass regions including the mass A ~ 40 region. 9 )~ n ) 
In the case of 44 Ti, the rotational band built upon the excited 0^ state at 1.905 MeV 12 ) 
was found to extend up to the J 7r =12 + state and was discussed to have superdeformed 
structure. 11 ) According to the shell model and mean-field-type calculations of 44 Ti, n )' 13 )" 15 ) 
this superdeformed structure is regarded as having a 4-particle-jump configuration, namely 
8p-4h configuration. The appearance of the superdeformed mean field in a very low excitation 
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energy region indicates the presence of the active dynamics of the mean-field formation. 

As is shown in the above discussions, 44 Ti is a typical nucleus which displays to us vividly 
the richness and profoundness of the nuclear many-body dynamics. Namely in 44 Ti we have 
the coexistence of the a-clustering dynamics and the mean-field dynamics, and also the 
coexistence of the normal deformed mean field and superdeformed mean field. Therefore it 
is an important and challenging task for the nuclear theory to describe properly and unifiedly 
these two kinds of coexistence phenomena of the nuclear structure dynamics displayed in 
44 Ti. We think that the proper description of these two kinds of coexistence phenomena 
in 44 Ti is indispensable for us in order to get the proper understanding of the neutron- 
proton correlation in the p /-shell region and also in the proton-rich region along N m Z 
line. In our recent study of 32 S, 16 ^' 17 ^ we have concluded that the superdeformed excited 
band which is predicted by many mean-field calculations 18 -* -20 -* to exist with its band-head 
+ state around 9 MeV above the ground state can be regarded also as being an 16 O + 16 O 
molecular band which is predicted to exist in the same excitation energy region from the 
i6q_i6q un jq Ue optical potential studies. 21 ^ 22 ^ This result of ours about 32 S demands us that 
we should study also in Ti the interplay and the relation between the cluster structure and 
the superdeformed mean-field structure. 

The purpose of this paper is to study the above-mentioned two kinds of coexistence 
phenomena in 44 Ti by using a single theoretical framework of the deformed-basis AMD 
(antisymmetrized molecular dynamics) . 2 )> 17 )> 23 ) The deformed-basis AMD is very suited for 
our present study because the AMD method is a kind of ab initio theory which enables us 
to study the clustering dynamics without assuming the existence of any clusters, 17 -* and also 
because, by the use of the deformed basis, the AMD can describe even the superdeformed 
mean field yielding almost the same quantitative results as the mean field theory. We will 
report that the AMD calculation by the use of the Gogny D1S force 24 -* surely reproduces 
and hence confirms the existence of the two kinds of coexistence in 44 Ti. Namely, firstly 
as for the confirmation of the coexistence of the a-clustering dynamics and the mean-field 
dynamics, we have obtained the K n =0~ band with prominent a+ 40 Ca structure and the 
K w =3~ band with prominent mean-field-like structure having a 1-particle-jump dominant 
configuration. Furthermore we have found the mixture of the a+ 40 Ca clustering character 
and the mean-field character in the K w =0 + ground band states. In the case of the ground 
state, the calculated a+ 40 Ca component is about 40%. When we judge from the usual 
expectation that in 44 Ti the mean field is formed more firmly than in 20 Ne and also the 
spatial-symmetry- breaking (hence a-cluster-breaking) spin-orbit force has stronger effects 
than in 20 Ne, it seems that we may doubt the a- clustering character in the ground band. 
However, the above-mentioned result of our calculation gives us an important theoretical 
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confirmation of the non-small a+ Ca clustering character of the ground band, which is 
consistent with the study by the use of the unique a- 40 Ca optical potential mentioned before. 
Just like we have obtained the K w =3~ band which is formed by the excitation of the single- 
particle degree of freedom from the mean-field-structure component contained in the ground 
band, we have obtained, on the other hand, an excited band which is formed by the excitation 
of the inter-cluster relative motion from the a+ 40 Ca-structure component contained in the 
ground band. This excited band with K n =0 + has been already reported to be observed 
experimentally and our AMD results are consistent to this experimental report. Secondly, as 
for the confirmation of the coexistence of the normal deformed mean field and superdeformed 
mean field, we have obtained, in addition to the normal deformation of the ground band and 
the K n =3~ band, the superdeformation of the K^—O^ band upon the 0^" state which has 4- 
particle-jump intrinsic configuration with large deformation parameter /3~0.5 and with large 
7 parameter 7~25°. Due to this triaxiality of the superdeformation we have also obtained a 
K w —2f superdeformed side band which we consider to correspond to the observed K n =2f 
band upon the 2jj~ state at 2.88 MeV. The triaxiality of the 44 Ti superdeformation was also 
reported in Ref. 15 - 1 

The contents of this paper are as follows: In the next section (section 2), we briefly 
explain the formalism of the deformed-basis AMD and the generator coordinate calculation in 
which the constrained deformed-basis AMD wave functions are superposed after the angular 
momentum projection. In section 3, we discuss the energy surfaces with good spins as 
functions of the deformation parameter (3 and the comparison of calculated energy spectra 
with experiments. The density distributions of the obtained intrinsic wave functions are 
discussed. The existence of the highly excited states with hyperdeformed structure is also 
suggested. In section 4, we discuss the obtained level scheme of 44 Ti which includes many 
kinds of the nuclear structure. Firstly, we explain the existence of the different kinds of the 
nuclear excitations in this nucleus together with making the comparison with the available 
experimental data such as the excitation energies and the E2 transition probabilities. Then 
we focus on the a+ 40 Ca clustering in the obtained level scheme. The a+ 40 Ca component in 
the ground band is analyzed and the excitations of the inter-cluster motion in the excited 
bands are discussed. Section 5 is for summarizing discussions. 

§2. Framework of the deformed-basis antisymmetrized molecular dynamics 

plus generator coordinate method 

In this section, the framework of the deformed-basis AMD+GCM is explained briefly. 
For more detailed explanation of the framework of the deformed-basis AMD, readers are 
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referred to references. 2 - ) ' 17 - ) The intrinsic wave function of the system with mass A is given 
by a Slater determinant of single-particle wave packets fi(r); 

@int = — F=det{</?i,</?2, ~;(Pa}, (2-1) 
^(r) = fa(r)x&, (2-2) 

where the single-particle wave packet (pi consists of the spatial fa, spin Xi an d isospin ^ 
parts. Deformed-basis AMD employs the triaxially deformed Gaussian centered at Z; as the 
spatial part of the single-particle wave packet. 



fa(r) oc expl - ^ v a (r a -Z ia ) 2 \, 

^ rr= r n y * 



a=x,y,z 

Xi = "iXt + PiXl, \(*i\ 2 + \Pi\ 2 = 1 

£i = proton or neutron. (2-3) 

Here, the complex number parameter Z, which represents the centroid of the Gaussian in 
the phase space takes independent value for each nucleon. The width parameters u x , v y and 
v z are real number parameters and take independent values for x, y and z directions, but are 
common to all nucleons. Spin part Xi is parametrized by aij and $ and isospin part £j is fixed 
to up (proton) or down (neutron). Zj, v x ,v y , v z and a,, f3i are the variational parameters and 
optimized by the method of frictional cooling. The advantage of the triaxially deformable 
single-particle wave packet is that it makes possible to describe the cluster-like structure 
and deformed mean-field structure within a single framework which is discussed in detail in 
reference. 2 - ) 

As the variational wave function, we employ the parity projected wave function as in the 
same way as many other AMD studies 

& = P^int = {l ^^nt, (24) 

here P x is the parity operator and 4>i nt is the intrinsic wave function given in Eq(2T). Parity 
projection makes possible to determine the different structure of the intrinsic state for the 
different parity states. 

Hamiltonian used in this study is as follows; 

H = f + V n + V c -f g , (2-5) 

where T and T g are the total kinetic energy and the energy of the center of mass motion, 
respectively. We have used the Gogny force with D1S parameter set as an effective nuclear 
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force V n . Coulomb force V c is approximated by the sum of seven Gaussians. The energy 
variation is made under the constraint on the nuclear quadrupole deformation by adding 
to H a constraint potential V cns t = v cns t{((3) 2 — /3g) 2 with a large positive value for v cns t- 
At the end of the variational calculation, the expectation value of V cnst should be zero in 
principle and in the practical calculation, we confirm it is less than O.lkeV. The optimized 
wave function is denoted by ^ ± ( / 9o). Here it should be noted that this constraint does not 
refer to the deformation parameter 7, which means that ^ ± ( / 9o) with positive [3q can be not 
only prolate but also oblate. 

From the optimized wave function, we project out the eigenstate of the total angular 
momentum J, 

*mk(0o) = Pmk^(Po) = P^P^intiPo). (2-6) 

Here Pmk * s the total angular momentum projector. The integrals over the three Euler 
angles included in the P(j K are evaluated by the numerical integration. 

Furthermore, we superpose the wave functions <£>j^ K which have the same parity and the 
angular momentum but have different value of deformation parameter /3q and K. Thus the 
final wave function of the system becomes as follows; 

<2>f = c n <^(/3 ) + (tfffM) + • • • , (2-7) 

where other quantum numbers except total angular momentum and parity are represented 
by n. The coefficients c n , c' n ,... are determined by the Hill- Wheeler equation, 

5(^1^1^) - E^l^)) = 0. (2-8) 

In order to analyze the obtained wave function ^MxiPo), we extract the single particle 
structure of <&i n t(f3o) by constructing the Hartree-Fock Hamiltonian from <Pi n t(/3o) as fol- 
lows. 25 ^ When the optimized wave function <£^ t (/3 ) = P ± -^ det{y?i, if 2, <Pa} is given, 
we calculate the orthonormalized basis <f) a which is a linear combination of the single-particle 
wave packets ifi, 

1 A 

(Pa = —= ^ (2-9) 

Here, \i a and Q a are the eigenvalue and eigenvector of the overlap matrix Bij = (ifi\ifj), 

A 

^ v PjjCja = l^aCia-, (2'10) 
J'=l 
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and it is clear that {<f) a } are orthonormalized from this relation. Using this basis set of 4> a , 
we calculate the Hartree-Fock single-particle Hamiltonian h a p which is defined as, 

A 

h a /3 = (</>a\i\<f>^) + ^(0q0 7 |t)|0/30 7 ) 
7=1 




(2-11) 



By the diagonalization of h a p, we obtain the single-particle energy e s and single-particle 
wave function <p s . 

A 

S ^h a pff3 S = e s f as , (2-12) 

(3=1 

A 
a=l 

We note that the single-particle energy e s and wave function <p s are obtained for occupied 
states but not for unoccupied states from this method. Furthermore, since the actual vari- 
ational calculation is made after the parity projection (the superposition of the two Slater 
determinants), it does not allow the naive interpretation of ± (/3o) by the single-particle 
picture. However, the single-particle structure obtained by this method is useful to estimate 
the particle- hole structure of the obtained wave function. From the parity nf = (<fi s \P ± \(() x ) 
and the the angular momentum (0 S |/ 2 |0 S ), we have estimated the particle-hole structure 
which is denoted as l nf\w structure' in the following. 

§3. Energy curves of the positive- and negative- parity states 

3.1. Positive parity states 

After the variation and angular momentum projection, we obtain the energy curves as 
functions of the matter quadrupole deformation parameter (3. Before discussing the low- 
lying states of 44 Ti, we study the features of the energy curves to clarify the characters 
of the wave functions on the energy curves. The positive-parity energy curve (figure 1) 
has three energy minima in each spin state and the wave functions at these minima have 
different particle-hole configurations. The wave functions around the minimum at /3=0.2 
have the Ohw configuration and they contribute to the ground band. Each spin state has a 
parity symmetric intrinsic density distribution and does not show the apparent existence of 
an a+ 40 Ca cluster correlation. The centroids of the single particle wave packets are gathered 
around the center of the nucleus, that is in contrast to the case of the 20 Ne ground band 
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in which single particle wave packets are separated into two parts (4+16) describing the 
spatially separated a and 16 O subunits. But as we will discuss in the next section, each spin 
states (except 10 + and 12 + ) around (3=0.20 contains a+ 40 Ca cluster correlation. Around 
/3=0.35, the change of the single particle configuration occurs and the wave function has a Afuv 
configuration in which four nucleons are excited into the p/-shell from 40 Ca core. This type of 
the configuration has its energy minimum around f3=0.5. Interestingly, the excitation energy 
of this Aftjjj state with respect to the Ohw state is quite small after the angular momentum 
projection. In the present result, it is only 1.7MeV in the case of the + state. The angular 
momentum projection lowers the energy of this state by about 6MeV. This state corresponds 
to the band head state of the superdeformed fT 7r =0^ band which starts from the 0^ state at 
1.9 MeV in experiments. Furthermore, this state has a triaxially deformed intrinsic density 
distribution (7=0.25) as is shown in figure 2. The triaxial deformation of this state have been 
also pointed out by the Hartree-Fock calculation 15 ) which is also free from any assumptions 
on the deformation symmetry of wave function. When we perform angular momentum 
projection, this triaxial deformation produces the K n =2 + band together with the K n =02 
band mentioned above. In figure 1, only the 2 + state is plotted among the K w =2 + band 
members to avoid the unsightliness of this figure. At larger deformation, the energy of the 
system increases rapidly, but again it has third local minimum around (3=0.7. Its energy is 
higher by about 19MeV than the Ohw state in the case of + state. Though the quadrupole 
deformation is not large enough, this state can be regarded as the hyperdeformed state, if 
we regarded the Ahw state around (3=0.5 as the superdeformed state. Because of the large 
parity asymmetry of the wave function, the system has the parity mixed single particle orbits 
and it makes hard to define the particle-hole structure of the usual single particle picture. 
It is notable that the density distribution of this state shows the possible relationship with 
the 16 0+ 28 Si molecular structure. Different from other states, the centroids of the single 
particle wave packets are separated into two parts (16 packets at left side and 28 packets at 
right side in figure 2 (e)) describing spatially separated 16 and 28 Si. 

3.2. Negative parity states 

Next, we discuss the negative-parity energy curves (figure 3) whose behavior is rather 
complicated than that of the positive-parity states. In the negative-parity states of this 
nucleus, the intrinsic wave functions have parity mixed single particle orbits. Therefore, 
we cannot conclude the definite t\w structure, though for the presentation we denote these 
structure with the ftw notation which is roughly estimated from the parity of the single 
particle orbit. In the small deformation region, we obtained the K w =3~ and K n =0~ bands 
and both have minima around (3=0.25. Roughly speaking, the intrinsic wave functions of 
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these states have lfuv structure. The K 7r =3~ band corresponds to the observed lowest 
negative-parity K n =3~ band, while the existence of the K n =0~ band in this energy region 
has not been established in experiments. In the medium deformed region around [3=0.4, we 
have obtained other K n =0~ and K n =3~ bands which are also not established in experiments. 
In these bands, 3hw and 5Hu) structures are mixed and they can be regarded as the lhw 
de-excited and excited states of the triaxial superdeformed state with the positive-parity, 
although their triaxiality are smaller than that of the positive-parity states. Their triaxial 
deformation are about 7=0.16° at most. 

Interestingly, we have also found the local minimum between the lfkv structure and 
3hu (5hcu) structure which is located around {3=0.27. This intrinsic state is obtained by 
performing the frictional cooling calculation using the a+ 40 Ca Brink wave function as the 
initial wave function. In this calculation, this structure continues to be stable during the 
cooling variational calculation after the parity projection to the negative-parity, although of 
course, the wave function considerably deviates from the initial Brink wave function through 
the cooling. It is contrasting to the case of the positive-parity state, in which the a particle 
is absorbed into the 40 Ca core even when we start the cooling from the a+ 40 Ca Brink 
wave function. This state corresponds to the theoretically suggested 5 - 1 ' 7 ) and experimentally 
observed 8 -* a+ 40 Ca cluster structure with the negative-parity. However, we will see in the 
next section that this band is strongly mixed with the above mentioned K w =0~ band {3fko 
and 5hu) after the GCM calculation and this mixing causes the fragmentation of the a+ 40 Ca 
rotational band states with negative-parity. In the largely deformed region, we found the 
counter part which originates from the parity asymmetry of the 16 0+ 28 Si structure obtained 
in the positive-parity state. The 16 0+ 28 Si structure with the negative parity gives rise to a 
K*=Q~ band. 

§4. Variety of the nuclear structure in the low-lying states of 44 Ti 

The calculated and observed level schemes of 44 Ti are presented in figure 5. In this 
section, we discuss the obtained level schemes together with other calculated quantities such 
as the E2 transition probabilities. Firstly, we discuss the different characteristics of the 
ground and excited states. Then we focus on the a+ 40 Ca clustering in this nucleus. We will 
see that the observed eight rotational bands shown in figure 5 are all reproduced well by the 
theory. 
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4.1. Low-lying level scheme 

4.1.1. K w =Of ground band 

The ground band member states mainly consist of the OHui structure and are obtained 
up to the 12 + state. The obtained binding energy of the ground state is -372.4MeV which 
underestimates the observed value by 3 MeV. The excitation energies of the ground band 
member states always underestimate the experimental values slightly. One of the reason of 
this deviation may be the pairing correlation which may not be included sufficiently in the 
present calculation, because it will act most strongly in the ground state. Of course, the 
insufficiency of the deformed-basis AMD wave function and/or the insufficiency of the Gogny 
D1S force to describe this band would be other reasons, since the excitation energy of the 
12 + state is too low to be explained by only the underestimation of the pairing correlation. 
Nevertheless we can appreciate the description of the ground band by the present model. 
The intra-band E2 transition probabilities of the ground band are reproduced well without 
any effective charge (tablel). It is to be noted that the calculated excitation energies show 
a similar characteristic behavior to the experimental values. Namely, in the low-spin states 
they have a rotational spectrum, but the 10 + and 12 + states show the diminish of the 
excitation energy. This change of the spectrum is closely related to the a+ 40 Ca clustering 
in the ground band and will be discussed together with its parity doublet partner, K^—O^ 
band in the next subsection. 

Table I. Observed (EXP) 26 )" 28 ) and calculated (AMD) intra-band E2 transition probabilities 
B(E2; J? -> Jj) [e 2 fm 4 ] within the ground band. For the comparison, o+ 40 Ca RGM (RGM) 7 ) 
and a+ 40 Ca(/ w ) OCM (OCM) 5 ) results are also shown. 



J i ~* J f 


EXP 


RGM 


OCM 


AMD 


2 + -0+ 
4 + -2+ 
6 + -4+ 
8 + -6+ 
10+ -> 8+ 
12+ -> 10+ 


120 ± 30 
280 ± 60 
160 ± 30 

14 > 
140 ± 30 
40 ±8 


107 

146 

140 

118 

75 

34 


166 
231 
212 
185 
138 
73 


142 

222 

167 

172 

99 

69 



4.1.2. K*=S[ and 0^ bands 

The Ifao excitation from the ground band produces the K w —0^ and K 7T =3± bands in 
the negative-parity. Though they have different K quanta, these two bands are obtained 
from the same intrinsic wave function. This induces the rather strong mixing between the 
K 7r =0 1 " band members and the natural parity band members of the band, and this 
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mixing enhances the even-odd staggering in the K w =3i band. The obtained K w =3i band 
corresponds to the experimental K n —3i band which is observed from the 3~ state at 3.17 
MeV to the 13~ state at 11.55 MeV. The lfiw nature of this band is consistent with the 
results of the shell model calculation. 11 ) We have assigned the calculated K n =0~[ band 
members to the observed 1~ (3.76 MeV), 3~ (3.94 MeV) and 5" (5.21 MeV or 5.31 MeV) 
states. In the present calculation, member states of this band are obtained up to 9~ state. 
Though this band is not sufficiently confirmed in experiments, the agreement between the 
calculation and the experiment seems to be satisfactory at least in the excitation energy. 

4.1.3. K^—O^ and 2 + superdeformed bands 

In contrast to the lftw single-particle excitation in the lowest negative-parity bands, the 
reconstruction of the mean-field occurs in the positive-parity. The first excited band in the 
positive-parity states is the K^—O^ band which mainly consists of the Ahv structure. As 
is discussed in the previous section, this 4fku structure has a triaxially deformed intrinsic 
wave function in which four nucleon are excited to the p/-shell from the 40 Ca core. Due to 
this many-particle excitation, the system has a quite different shape of the mean-field and 
is stabilized by the formation of the elongated superdeformed structure. After the angular 
momentum projection and the GCM calculation, the excitation energy of this band becomes 
quite small. It starts from the 0^ state at 1.75 MeV and reaches the 10 + state at 12.1 
MeV in the present calculation. The obtained excitation energies show good agreement with 
experiments in low-spin states. As the angular momentum increases, the observed spectrum 
deviates from the rotational one, while the calculated spectrum shows the rotational charac- 
ter. This may be because of the absence of the 2ftw structure in our calculation. In the shell 
model calculation, 11 ^ it is shown that the 2htv structure becomes dominant as the angular 
momentum increases in this band and the excitation energies of 10^" and 12J, states are 
much lower than those expected from the rotational spectrum. In our calculation, we have 
superposed the wave functions on the energy curve and the 2ftw structure has not appeared 
on the energy curve. Of course, the obtained wave function in the present calculation is not 
of the pure Ohu and 4hw structure, but also contains 2hw structure to some extent. However 
the description of the 2hu structure seems to be not sufficient in the present calculation. 
This may also explain the absence of the 12^" state in our calculation. 

Since the 4ftw structure has a triaxially deformed intrinsic wave function, K n =2 + band 
adjoins K^—O^ band. The member states of this band are obtained from the 2g state (3.2 
MeV) up to 10 + (11.3 MeV) state. In the experiment, the possible existence of the low-lying 
K n =2 + band which starts from the 2 + state at 2.89 MeV has been long known, and more 
recently the 8 + , 10 + and 12 + states which possibly belong to this band have been observed 
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Table II. Observed (EXP) 26 )~ 28 ) and calculated (AMD) intra- and inter-band E2 transition prob- 
abilities B(E2;J[ -> JJ) [e 2 fm 4 ] within the K n =0^ and K 7T =2+ bands. For the comparison, 
a+ 40 Ca(/ 7r ) coupled-channel OCM (OCM) 5 ) results are also shown. 





EXP 


OCM 


AMD 


^2 ~^ U 2 


oon 4- c;n 

ZZU It ou 

268 ± 50 


1 %l7 
10 1 

268 


ozu 
361 


7^=2+^7^=2+ 


EXP 


OCM 


AMD 


3j — > 2+ 
4+^3+ 


< 590 
1 7C+100 

1 '°-60 

< 785 ± 650 


185 
148 
11 


298 
220 
302 




EXP 


OCM 


AMD 


2+^0+ 


43 < 


3.04 


24 



together with their 7 transitions to the K n —02 band, though the spin-parity of the (3 + ) state 
is not fixed and the 5 + , 6 + , 7 + , 9 + and 11 + states have not been observed yet. The obtained 
K 7T =2 + band corresponds to this K w =2f side band. The calculated inter- and intra- B(E2) 
values within K^—O^ and K n =2f bands (table II) are larger than the experimental values, 
but show similar trend qualitatively. 

We note that these two band were also obtained up to 4 + states by a semi-microscopic 
cluster model 5 -* assuming the a+ 40 Ca(0 + ) and a+ 40 Ca(2+) cluster structure. This result also 
supports the 4fuv structure of these bands and means the same origin of the 7C"=0 + and 
K 7T =2 + bands, since both 40 Ca(0+) and 40 Ca(2^) states are regarded as having the ct+ 36 Ar 
cluster structure in their calculation. We also note that their results implies that the K^—O^ 
and K 7r —2 + bands are possibly regarded as having the 2a+ 36 Ar cluster structure. However, 
in the case of our AMD result the relation between the superdeformed state and the 2a+ 36 Ar 
structure is unclear. Indeed, in our calculation this state has a quite large contribution from 
the spin-orbit force (about -27 MeV in the case of the + state) which is much larger than 
that expected in the case of the 2a+ 36 Ar structure (about -14 MeV in our calculation). 

4.1.4. Hyper deformation and 16 0+ 28 Si molecular structure 

About 19 MeV above the ground state, we have obtained the K w =0 + band which will 
corresponds to the f 8 g 4 state obtained in the Hartree-Fock calculation 15 ) with large defor- 
mation /9~0.7 and it can be classified as a hyperdeformed state. However, as was shown 
in the previous section, this state has an 16 0+ 28 Si intrinsic structure and it implies that 
it will be more suitable to understand this state as the 16 0+ 28 Si molecular structure than 
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to regard as the hyperdeformed state. Indeed this band has its counterpart in the negative 
parity, K n =0~ band which also has an 16 0+ 28 Si intrinsic wave function. It may be possible 
to regard these two bands as the parity doublet bands. These bands can have something to 
do with the possible existence of the 16 0+ 28 Si resonant states. 

4.2. a+ 40 Ca cluster structure 

In the previous subsection, we have shown that many kinds of the nuclear excitation 
appear in the 44 Ti and their description by the present model shows satisfactory agreement 
with experiments. Besides these structures, the existence of the a+ 40 Ca cluster structure in 
44 Ti has long been discussed by many authors. In particular, the semi-microscopic cluster 
model which uses the a+ 40 Ca optical potential and the a+ 40 Ca resonating group method 
calculations 3 )' 4 )' 7 ) have shown the a+ 40 Ca nature of the ground band and the existence of 
the excited N=13, 14 and 15 bands of the a+ 40 Ca structure negative-parity. Here N denotes 
the principal quantum number of the relative motion between a and 40 Ca, and is defined 
as N=2n + L with n and L denoting the number of nodes and the angular momentum of 
the a-Ca relative wave function, respectively. N=12 is the lowest Pauli allowed state which 
corresponds to the ground state. In this section, we discuss the a+ 40 Ca clustering in the 
obtained states on the basis of the microscopic model, deformed-basis AMD, which can also 
reproduce the excitation mechanism other than the clustering. To clarify our argument, we 
show the results of two GCM calculations which are obtained by superposing the Ohio wave 
functions and a+ 40 Ca wave functions on the energy curves in Fig. 1 and 3 (Qftw+(a+ m Ca) 
GCM) and by superposing all of the obtained wave functions (full GCM). The result of the 
0fuv+(a+' l0 Ca) GCM shows the existence of the non-small a+ 40 Ca component in the ground 
band and the lhcu, 2f\w and 3ftw excitation of the relative motion between a and 40 Ca in the 
N=13, 14 and 15 bands. The comparison with experiments should be made by the full GCM 
calculation. We will see that, due to the strong mixing with the 3ftw structure which does 
not have the a+ 40 Ca component, N=13 and 15 bands are fragmented into several states in 
the full GCM calculation. 

4.2.1. a+ 40 Ca component in the ground band 

In the present calculation, the ground band is the only candidate of the lowest a+ 40 Ca 
band among the obtained K 7T =0 + bands, since other obtained K 7T =0 + bands have 2hu or 
ihuj excited structure of the 40 Ca core. We found that the ground band member states have 
non-small amounts of the a+ 40 Ca component except 10| and 12+ states. The method of the 
evaluation of the a+ 40 Ca component in the AMD wave function is explained in the appendix. 
In the case of the ground state obtained by the full GCM calculation, it amounts to about 
39% (table III). However, it is not so large as in the case of the a+ 16 component in the 
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ground state of 20 Ne, 2 - ) in which it amounts to about 70%. As is stated in the introduction, 
due to the increasing core mass and the strong effect of the spin-orbit force, we cannot 
expect the prominent a+ 40 Ca cluster structure in the ground band, which is different from 
the cases of the sci-shell nuclei such as 16 and 20 Ne. Indeed, the spin-orbit force has a non- 
small expectation value (-9.5 MeV in the case of the ground state) which must be zero in the 
case of the a+ 40 Ca structure. Therefore, we can conclude that the a+ 40 Ca cluster structure 
is considerably distorted by the formation of the deformed mean-field and the strong effect 
of the spin-orbit force. 

To understand the distortion of the a+ 40 Ca cluster structure in the ground band, we 
first decompose the wave function of the ground state obtained by the full GCM calcula- 
tion into the a+ 40 Ca cluster component and other components orthogonal to it. Then we 
decompose the a+ 40 Ca cluster component into the SU3 irreducible representations. The 
definitions and explanations of this decomposition are presented in the appendix. Similarly 
to the ground state of the 20 Ne, the ground state of 44 Ti contains the highly excited (N, 0) 
components (figure 7). For example, (N m i n +8, 0) = (20, 0) component which corresponds to 
the 8hw excitation in the spherical harmonic oscillator basis still amounts to about 3%. This 
coherent contribution from (N, 0) states is specific to the cluster correlation. However, we 
must note that the magnitude of these contribution is quite small compared to the 20 Ne. 
In particular, the contribution from the lowest Pauli-allowed state, ((12,0) state in the case 
of the 44 Ti and (8, 0) state in the case of the 20 Ne) is greatly diminished in the 44 Ti. This 
means that the cluster structure is considerably dissolved inside the nucleus, since smaller 
N means the smaller inter-cluster distance. And the residual coherent contribution from 
higher (TV, 0) components means that the a+ 40 Ca clustering survives in the surface region, 
but its contribution is not so large as in 20 Ne. In other words, the binding mechanism of the 
ground band of 44 Ti may be rather dominated by the mean-field picture and the a+ 40 Ca 
cluster correlation shows up more clearly in the surface region of the nucleus. It is to be 
noted that the ground band still has a considerable amount of the a spectroscopic factor S a . 
Since the smaller (N, 0) components less contributes to the S a due to the smallness of the fi^ 
(see Eq.(A-9) in the appendix), the dissolution of the a+ 40 Ca component in the smaller dis- 
tance (smaller (N, 0) components) does not affect so much to the S a . The dominance of the 
mean-field-like component becomes more prominent as the angular momentum increases. In 
particular, the 10^~ and 12^ states has quite small amount of a+ 40 Ca components. Instead, 
in these states, the intrinsic wave functions has an oblately deformed shape. This structure 
change along the yrast line causes the irregular behavior of the excitation energies of lOf 
and 12 1" states. The change of the deformation to the oblate shape near the band terminal 
is due to the nucleon spin alignment to the rotation axis and was also observed in the AMD 
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calculation of the 20 Ne ground band. 29 -* 

4.2.2. K n =% band with a+ 40 Ca structure 

It must be noted that the considerable distortion of the a+ 40 Ca cluster structure in the 
ground band does not mean the less importance of the a+ 40 Ca clustering in this nucleus. 
In the previous section, we have shown that the lfuv single-particle excitation from the 
ground state produces the K w —3i and K n =0^ bands. From these state, the single-particle 
excitation is regarded as the one of the basic degrees of the freedom of the excitation in this 
nucleus. In the same sense, the a+ 40 Ca clustering is also regarded as one of the basic degrees 
of freedom of the excitation, since the K n —0~ band in which the system has a considerable 
amount of the a+ 40 Ca component is obtained in the present microscopic model calculation. 
As is shown in the previous section, the wave function on the energy curve which forms 
the local minimum around (3=0.35 shows a existence of the a+ 40 Ca clustering in this state. 
Indeed, when we perform the 0hu+(a+' l0 Ca) GCM calculation, it produces the K w =0~ band 

Table III. Amount of the o+ 40 Ca component Wj and a spectroscopic factor S a of the obtained 
states (see text and appendix). The full GCM results of N=13 and 15 band members and the 
experimental value of N=13 band members show the sum of the fragmented states. 
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(N=13 band) which starts from the 1~ state at 8.9 MeV and the member states of this band 
have considerable amount of the a+ 40 Ca component (tablelll). In the case of the 1~ state, 
it amounts to about 60%. This band lies in the same energy region as the experimentally 
observed states. The decomposition of the relative wave function between a and 40 Ca shows 
the growth of the a+ 40 Ca clustering in this band. Compared to the ground band, the amount 
of the lowest component ((12,0) in the case of the ground band and (13,0) in the case of 
the K n =02 band) considerably increases together with the coherent contribution from the 
higher (N, 0) states and they lead to the increase of the S a in this band (table III). We can 
say that since the ground state contains the non-small amount of the a+ 40 Ca component 
by about 40%, the nucleus exploits this degrees of the freedom in the nuclear excitation. 
However, when we perform the full GCM calculation, the member states of this band are 
fragmented into several states because of rather strong mixing with 3fiuj structure which 
does not have a+ 40 Ca component. In figure 6, the excitation energies of these fragmented 
states and their averaged energies are plotted. The present model space is not large enough 
to be directly compared with the observed fragmented states. But, at least, we may be able 
to say that one of the origin of the observed fragmentation is the mixing of the a+ 40 Ca 
structure with the 3hv non a-clustering structure. 

4.2.3. Excited K 7r =0 ± bands with a+ 40 Ca structure 

The a+ 40 Ca clustering as a degree of freedom of the nuclear excitation is reinforced by 
the existence of the excited state in which the relative inter-cluster motion is excited and 
have the principal quantum numbers iV=14 and 15. When we perform the 0fkJ+(a+ 40 C&) 
GCM, we obtain additional two rotational bands in positive- and negative-parities. These 
two bands starts from 0j state at 15.1MeV and 1 3 state at 16.9 MeV, respectively, and 
they contain considerable amount of the a+ 40 Ca component. We confirmed that these two 
bands corresponds to iV=14 and 15 bands by investigating the relative motion (figure 7). 
Namely, iV=14 band has a dominant (14,0) component and iV=15 band has a dominant 
(15,0) component. The iV = 15 band has a largest amount of the a+ 40 Ca component among 
7V=12,13,14 and 15 bands. This is because of the largest inter-cluster distance of this band 
which is due to the Afkj excitation of the relative motion. 

When we perform the full GCM calculation, the iV = 15 band members also couple with 
3fru structure and are fragmented. But the energy centroids (averaged energy) are not so dif- 
ferent as the excitation energies of the non- fragmented states obtained by the Oftu;-|-(a:-|- 40 Ca) 
GCM. The iV=14 band has not been fragmented in the present calculation, though its frag- 
mentation is observed in the experiment. We consider that this is due to the absence of 
the 2fkv structure in our calculation. As is already mentioned, the 2fkv structure does not 
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appear on the obtained energy curve. And we guess that 2ftw structure is most likely to be 
coupled with the A=14 band, since their principal quantum numbers are equal. 

§5. Summary 

We have studied the low-lying level structure of 44 Ti by using the framework of the 
deformed-basis AMD with the Gogny D1S force. We have shown that many kinds of mutually 
very different nuclear excitation appearing in this nucleus are very well reproduced unifiedly 
in a single AMD framework. Namely, the single-particle excitation in the K n =3i and 
K n —0i bands, the formation of the triaxial superdeformed states in the K v =^2 and K n =2 + 
bands, and the a+ 40 Ca clustering in some K n =0 ± bands. The possible existence of the 
16 0+ 28 Si molecule-like structure is also suggested in the present calculation. The triaxial 
superdeformed state suggested by some Hartree-Fock calculations are confirmed and it has 
been shown that this superdeformed state explains the low-lying K^—O^ and 2 + bands, 
which suggests the interpretation of the low-lying states of these bands as the triaxial rotor. 
Among these state, we have focused on the a+ 40 Ca clustering in the low-energy region. We 
found that the a+ 40 Ca cluster component in the ground band is considerably dissolved by 
the formation of the deformed mean-field and the spin-orbit force. But, since the ground 
state still contains non-small amount of the a+ 40 Ca component by about 40%, it is possible 
to excite the nucleus by using the clustering as a degree of the freedom of nuclear excitation. 
This excitation leads to the existence of the i^ 7r =0 ± bands in which the relative motion 
between a and 40 Ca have the principal quantum number N=13, 14 and 15. However, due to 
the coupling to 3fkv structures, N=13 and 15 bands members are fragmented into several 
states. When we average their energies for each spin, these bands lie in the same energy 
region of the observed a+ 40 Ca bands with N=13, 14 and 15. 
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Appendix 



The definition and the explanation of the a+ 40 Ca component in the AMD wave function 
and its decomposition into the SU3 irreducible representations are presented. The reader is 
directed to references 30 )' 31 ) for more details. 

A.l. a+ 40 Ca component in the AMD wave function 

The a+ 40 Ca system is generally expressed by the RGM-type wave function. 

<-4o Ca = ^{ Xj (OX/o(r)0(a)0( 4O Ca)}, (A-l) 
n = ^40! • 4!/44!, tt = (-) J . (A-2) 

Here, A is the antisymmetrizer, r is the relative coordinate between a and 40 Ca, 0(a) and 
0( 4O Ca) are the normalized internal wave function of the clusters. Xj( r ) is the radial wave 
function of the relative motion between a and 40 Ca, and is so normalized that &a+ 40 c& is 
normalized to unity. 

The normalized deformed-base AMD+GCM wave function <P{^ of 40 Ca is divided into 
the a+ 40 Ca component and the residual part <P R . 

< = ^f +4 o Ca + VT=tf*£ . ( A3) 

is also normalized and orthogonal to the ^+ 40 c a - We introduce the projection operator 
Vl which projects out the a+ 40 Ca component from the $m , 

V L $ti = a^ 40Ca 

= an o ^{xj(r)F JO (r)0(«)0( 4O Ca)}. (A-4) 

The practical formula of Vl used in this study is given later. Using this projection operator, 
the squared amplitude of the a+ 40 Ca component (amount of the a+ 40 Ca component) Wj 
of $m is written as, 

Wj = \a\ 2 = (^\V L \<Pij). (A5) 

A. 2. Decomposition of the a-/- 40 Ca component into SU3 irreducible representations 

If the internal wave functions 0(a) and 0( 4O Ca) are expressed by the H.O. wave functions 
with oscillator parameter v, Eq. (AT) and (A-4) are easily decomposed into SU3 irreducible 
representations. When we expand Xj( r ) by the radial wave function of the H.O. R^j^r, 7) 
with N=2n + J and the width parameter 7 = as 

Xj(r) = e NjRNj(r, 7), (A-6) 

N 
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the norm of Eq. (AT) is written as 



N 
_40, 



and therefore, the squared amplitude of the a+ Ca is rewritten as 

WT = \a\ 2 = \a\ 2 J2^Ne 2 NJ . (A-8) 



AT 

and the a spectroscopic factor S a is defined as, 



Si^\a\ 2 J2^e 2 Nj- (A-9) 



Here, the eigenvalue of the RGM norm kernel /xjv is defined as 

fi N =( J R 7VJ (r)F J o(f)0(a)0( 4O Ca)| 

^{i^(r)y JO (r)0(a)0( 4o Ca)}>, (A-10) 

and it has a non-zero value for N > 12 in the case of the a+ 40 Ca system. The Eq. (A-10) 
means that the quantity \a\ 2 ii^e 2 NJ is equal to the amount of the SU^(N,0) irreducible 
representations which is contained in ct^+ 40 Ca> srnce both <j)(a) and 0( 4o Ca) are SU3 scalar 
and therefore, the a+ 40 Ca system which has a relative motion xj( r )=RNj(r, 7) belongs to 
SU 3 (N,0). 

A. 3. Approximation of the projection operator Vl 

To evaluate Wj and ejvj, we have approximated the projection operator Vl by the set 
of the orthonormalized angular-momentum-projected Brink-type wave functions. We start 
from the Brink-type wave function (p B (Ri) in which a and 40 Ca are located at the points 
§Ri) and (0,0,-^ 



(0, 0, gi^) and (0, 0, -A-Ri), respectively, 



MRi) =n ^l{^(a,^R l ),^( 40 Ca,-^R i )}, (ATI) 
^ = (0,0,^). " (AT2) 

For simplicity, tp(cx) and (/?( 40 Ca) are described by the AMD wave functions which are ob- 
tained by the variational calculation without any constraint potential and with the fixed 
spherical width parameter v— 0.150. And we have assumed that these wave functions ap- 
proximate well the SU 3 limits. By separating the center-of-mass wave function uj(Xq), the 
internal wave function ip c (Ri) is written as 

<p B (Ri) =cu(X G )<pc(Ri), (A-13) 
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MRi) = n A{r(r, K t , 7 )0(a)0( 4O Ca)}, (A-14) 

f(r,R„ 7 ) = ^) 3 e-^' 2 , (A-15) 
/44-2z/\ 3/4 

u(Xa) = \^ L ^) exp(-44i/X£), (A-16) 



44 

I 



44^ " ' 44 

i=i 

where -T(r, Ri, 7) is the relative wave function between a and 40 Ca of the Brink-type wave 
function. The angular-momentum-projected Brink-type wave function is obtained 

from the <f>c(Ri), 

<p J c (R i )=q(P^ c (Ri)- (A-18) 
P ^ is the angular momentum projector and qf is the normalization factor. The orthonor- 
malized set of the angular- momentum-projected Brink-type wave functions Cp J a is described 
by the linear combination of the </?^(i?j), 

and the coefficients p a and Cj a are the eigenvalues and eigenvectors of the overlap matrix 
fly H^MbcTO), 

^ ] Bjj c ja — Pa c ia- (A-20) 
3 



If the number of the basis wave functions (p a is taken enough, the projection operator Vl is 



xJ 

1UUV.UU110 ^ 

approximated by the set of 

= E S 7*VcW><Vc(^)l- (A-21) 

u 

In the present calculation, we have employed 19 Brink-wave functions in which Ri is taken 
as = 1.0, 1.5, 2.0,. ..,10.0 fm. 

Using Eq. (A-21), Eq. (A-4) and Eq. (A-5) are rewritten as 

Vl&jS = an o ^{xj(r)F7o(r)0(a)0( 4O Ca)}, 

^E^^c^OlObcW), (A-22) 
y 



x«>£(^)>- (A-23) 
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The expansion cofficient ejy of Xj{ r ) is obtained from Eq. (A-22). We expand wave 
function of the relative motion 7(1*, Rj, 7) of Eq. (A-24), 



r(r,R i ,7) = e-i^^ 

NJm ^ 

X 



]/ ^^ A J Y Jm(Ri)RNj(r, 7 )V>m(r) }, (A-24) 
.Af _ / \(N-J)/2 (2J+1)-N\ 

Aj " (j y (AT — J)!! • (N + J + 1)!! ' (A25) 
And using Eq. (A- 14) and (A-24), Eq. (A- 18) is rewritten as 

V&Ri) = n o ^l{x5 ) (r)F7o(r)0(a)0( 4O Ca)}, (A-26) 

Xf(r) = £ ^^i2 w (r). (A-27) 

N ^ iv. 

By inserting Eq. (A-26) into Eq. (A-22), we obtain Xj{ r ) as a superposition of Xj\ r )i 

Xj(r) ^^^{R^Bjlxfir). (A-28) 

ij 

Inserting (A-27) into Eq. (A-28), is given as 

^ = \a n j Y,\^ , c(R 3 )Wm)bj 1 1 



a 



N 



x?/e -i«f(2jl|. (A . 29) 
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Fig. 1. Energy curves for each spin state as functions of the quadrupole deformation parameter 
p. + , 4 + , 8 + and 12 + states are shown for the presentation. The 2 + state which is the band 
head state of K w =2 + band is shown by the dotted line. The parity-projected intrinsic state 
(before the angular momentum projection) is also shown by the dashed line. 
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Fig. 2. Density distributions of the intrinsic wave functions obtained by the variational calculation 
after the projection to positive parity, (a) shows the density distribution for the OHuj configu- 
ration, (b), (c) and (d) show that for the Ahw configuration plotted from different angles, (e) 
shows that for the 16 0+ 28 Si configuration. Black points in these figures show the forty-four 
centroids of the single particle wave packets (5ft(Zj)) in the coordinate space. 
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Fig. 3. Energy curves for each spin state as functions of the quadrupole deformation parameter (3. 
Only the band head states are shown for the presentation. The parity-projected intrinsic state 
(before the angular momentum) is also shown by the dashed line. 
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Fig. 4. Density distributions of the intrinsic wave functions obtained by the variational calculation 
after the projection to negative parity, (a) and (c) shows the density distribution for the lhw 
and 3Hlj configurations, respectively, (b) shows that for the a+ 40 Ca configuration, (d) shows 
that for the 16 0+ 28 Si configuration. Black points in these figures show the forty-four centroids 
of the single particle wave packets (3ft(Zj)) in the coordinate space. 
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Fig. 5. Observed and calculated level scheme of 44 Ti. Only the observed states which correspond 
to the calculated states are shown. In the calculated results, the main component of each band 
is shown below the band head state (ex. Ohuj) as a guide. The states plotted by the bold lines 
in both experiment and theory correspond to the a+ 40 Ca cluster states. These states except 
of N=14 band in the theoretical result are shown by the averaged energy of the fragmented 
states (see next subsection). 
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Fig. 6. Observed 8 ) and calculated (0/kJ+(a+ 40 Ca) GCM and full GCM ) states of 44 Ti which have 
the a+ 40 Ca cluster structure. The bold lines indicates the energy centroid of the experimental 
(A/=13,14 and 15) and theoretical (full GCM, N=13 and 15) fragmented states. The dotted 
lines indicate the energy of the fragmented states obtained by the full GCM calculation. 
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Fig. 7. The decomposition of the a+ Ca component into the SU3 irreducible representations. 
These results are calculated from the Ofou+(a+ 40 Ca) GCM wave functions. For comparison, 
the result of the ground band of 20 Ne is also shown. N m i n indicates the lowest Pauli-allowed 
states. N m i n =12 and N m i n =8 for 44 Ti and 20 Ne, respectively. 
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